Observation of Dirac Holes and Electrons in a Topological Insulator 
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We show that in the new topological-insulator compound Bi1.5Sbo.5Te1.7Se1. 3 one can achieve a 
surfaced-dominated transport where the surface channel contributes up to 70% of the total conduc- 
tance. Furthermore, it was found that in this material the transport properties sharply reflect the 
time dependence of the surface chemical potential, presenting a sign change in the Hall coefficient 
with time. We demonstrate that such an evolution makes us observe both Dirac holes and electrons 
on the surface, which allows us to reconstruct the surface band dispersion across the Dirac point. 

PACS numbers: 73.25.+i, 71.18.+y, 73.20.At, 72.20.My 



The three-dimensional (3D) topological insulator (TI) 
hosts a metallic surface state that emerges due to a non- 
trivial Z2 topology of the bulk valence band [l], Q . This 
peculiar surface state offers a new playground for study- 
ing the physics of quasiparticles with unusual dispersions, 
such as Dirac or Majorana fermions 0,0]. However, most 
of the known TI materials are poor insulators in their 
bulk, hindering transport studies of the expected novel 
surface properties 0, Last year, we discovered that 
the TI material Bi2Te2Se_ (BTS) presents a high resis- 
tivity exceeding 1 51cm [5(, which made it possible to 
clarify both the surface and bulk transport channels in 
detail. Also, we found that in our BTS sample the surface 
channel accounts for ^6% of the total conductance. For 
this BTS compounds, Xiong et al. recently reported a 
higher resistivity in the range of 5-6 ficm, together with 
pronounced surface quantum oscillations which possibly 
signify fractional- filling states . Since the source of the 
residual bulk carriers in BTS is the acceptor states 
reducing the number of anti-site defects working as accep- 
tors in this promising material is an important challenge 
for the advancement of the TI research. 

In this work, we tried to optimize the composition of 
BTS by reducing the Te/Se ratio and introducing some 
Sb into Bi positions Q, while keeping the ordering of 
the chalcogen layers as in BTS (Fig. 1(a), Q). The 
X-ray powder diffraction patterns shown in Fig. 1(b) 
demonstrate that the chalcogen ordering is still present 
in Bi1.5Sbo.5Te1.7Se1. 3 (BSTS), and we focus on this com- 
pound in this Letter. We found that in BSTS one can 
achieve an even larger surface contribution (up to 70%) 
than in BTS. We also found that the surface state of 
BSTS changes with time, and, intriguingly, we observed 
quantum oscillations from Dirac holes, the negative en- 
ergy state of the Dirac fermions, as well as those from 
Dirac electrons in the same sample at different time 
points. We show that this time evolution can be used to 
reconstruct the surface band structure across the Dirac 
point, providing a unique way to determine the disper- 
sion relation of the surface state. 

Single crystals of BSTS were grown by melting high 
purity (6N) elements of Bi, Sb, Te, and Se with a mo- 
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FIG. 1: (Color online) (a) Basic structure unit of 
Bii.gSbo.sTei.TSer.a (BSTS). (b) Comparison of the X-ray 
powder diffraction patterns of BSTS, Bi2Te2Se, and Bi2Te3. 
Dashed lines indicate the positions of the peaks characteristic 
of the ordering of Se and Te (Te/Se) layers, (c) Temperature 
dependence of p xx measured repeatedly in time in a cleaved 
30-/iin-thick BSTS sample, (d) Temperature dependence of 
the low-field Rh , presenting a sign change with time. Dashed 
line represents the Arrhenius-law fitting to the data above 150 
K, which is also shown in the inset. 



lar ratio of 1.5:0.5:1.7:1.3 at 850°C for 48 h in evacuated 
quartz tubes, followed by cooling to room temperature 
over one week. For transport measurements, crystals 
were cut along the principal axes, and cleaved just before 
the measurements. Ohmic contacts were prepared by us- 
ing room-temperature cured silver paste. The resistivity 
p xx and the Hall resistivity p yx were measured simulta- 
neously by a standard six-probe method Q by sweeping 
the magnetic field B between ±14 T at 0.3 T/min while 
stabilizing the temperature T to within ±5 mK. 

Freshly cleaved thin samples were used for studying 
the surface transport in BSTS. As shown in Fig. 1(c), 
p xx in BSTS sharply increases upon lowering tempera- 
ture from 300 K, which is characteristic of an insulator, 
but it saturates below ~100 K due to the metallic sur- 
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FIG. 2: (Color online) (a,b) dp xx /dB vs 1/B± measured in 
tilted magnetic fields in the p state (5 h) and the n state (730 
h). All curves are shifted for clarity. Insets show the Fourier 
transforms of the data at (9=0°. (c,d) Temperature depen- 
dences of SdH amplitudes for 9=0° shown in (a) and (b) and 
their theoretical fittings. Insets show AR XX vs. 1/B after 
subtracting a smooth background, (e) Landau-level fan dia- 
gram obtained from the oscillations in AR XX measured at T 
= 1.5 K and 8 — 0°; minima and maxima in AR XX correspond 
to riosc and n osc + o i respectively. 



face transport as well as the bulk impurity-band (IB) 
transport [5j|. This behavior is essentially the same as in 
BTS. What is peculiar in BSTS is that p xx {T) increases 
slowly with time and, furthermore, the low-temperature 
Hall coefficient Rh changes sign with time in thin sam- 
ples. As an example, Figs. 1(c) and (d) show p xx {T) 
and Rh(T) data of a 30-/im-thick BSTS sample, mea- 
sured repeatedly in about one month. In contrast to the 
relatively small change in p xx [Fig. 1(c)], Rh at low- 
temperature exhibits rather drastic evolution [Fig. 1(d)] 
from positive to negative, whereas Rh at high temper- 
ature is essentially unchanged with time. This suggests 
that the source of the time dependence resides in the 
surface channel. In passing, Rh{T) above 150 K is pos- 
itive and demonstrates an activated behavior [shown by 
the dashed line in Fig. 1(d)], signifying the thermal ex- 
citation of holes into the valence band with an effective 
activation energy of about 60 meV. This is similar to 
what we observed in BTS [5j]. 

To understand the nature of the time evolution, the 
Shubnikov-de Haas (SdH) oscillations were measured in 
the aforementioned 30-/im-thick sample twice, 5 h after 



cleavage (called p state) and 725 h later (called n state), 
between which Rh changed sign (the sample was kept 
at 300 K in air). In BSTS, the oscillations do not fade 
out even after long exposure to ambient atmosphere, as 
opposed to other TI materials like Sb-doped Bi 2 Se3 [lOj . 
The SdH oscillations were clearly observed in p xx (B) in 
our samples, but they were not clear in p yx (B), so the 
following SdH analysis is restricted to p xx {B). Figures 
2(a) and (b) show dp xx /dB for both the p and n states, 
plotted as a function of 1/B± (= 1/B cos 9), where 9 is 
the angle between B and the C3 axis. Several equidis- 
tant maxima and minima are easily recognized, and im- 
portantly, the positions of maxima and minima depend 
solely on B±, which signifies a 2D character of the ob- 
served oscillations. Insets show the Fourier transform 
of the oscillations taken at 9=0°. Two frequencies are 
clearly seen in the p state, but the second one (60 T) is 
probably a harmonic of the primary frequency F = 30 T. 
On the other hand, only the primary F = 50 T is seen in 
the n state. The averaged Fermi wave number kp is ob- 
tained by using the Onsager relation F = (hc/2ire)irk F , 
resulting in kp = 3.0xl0 6 cm -1 and 3.9xl0 6 cm -1 for 
the p and n states, respectively. This corresponds to the 
surface hole (electron) concentration of 7.2xlO n cm~ 2 
(1.2xl0 12 cm -2 ) for a spin-filtered surface state. Fitting 



the standard Lifshitz-Kosevich theory [ll| to the tem- 
perature dependences of the SdH amplitudes [Figs. 2(c) 
and (d)] gives the cyclotron mass m c of (0.10±0.01)m e 
for holes and (0.075±0.003)m e for electrons (m e is the 
free electron mass). Also, from the independence of the 
SdH amplitudes [insets of Figs. 2(c) and (d)] one can 
obtain the scattering time r of 5.8xl0~ 14 s (4.2xl0~ 14 
s) for holes (electrons) through the Dingle analysis. 

From the measured values of kp and m c , one obtains 
the effective Fermi velocity v F (= hkp/m c ) of 3.5xl0 5 
m/s and 6.0xl0 5 m/s for holes and electrons, respec- 
tively. Now we discuss that this difference between the 
two v F bears crucial information regarding the Dirac 
cone: If the surface state consists of ideal Dirac fermions, 
the Fermi velocity is independent of k and is particle- 
hole symmetric. However, the energy dispersions of the 
surface states in Bi-based TI materials generally deviate 
from the ideal Dirac cone, and it can be described as [12j 
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where Edp is the energy at the Dirac point (DP), vf is 
the Fermi velocity at the DP, and m* is an effective mass 
that accounts for the curvature in E(k). The effective 
Fermi velocity v* F reflects the local curvature in E(k) and 
can be expressed as v* F {k) = (dE/dk)/h = Vp + hk/m*. 
The p and n states correspond to the situations where the 
Fermi energy Ep is below or above the DP, respectively, 
and the time evolution of Rh is a manifestation of the 
time-dependent change of the surface chemical potential. 
By using the kp and v* F values obtained for the p and 
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n states, we can solve simultaneous equations to obtain 
vp = 4.6xl0 5 m/s and m* = 0.32m e . This mean that 
the time evolution of the transport properties allowed us 
to do a "spectroscopy" of the surface state to determine 
its dispersion, from which we can further estimate the 
position of Ep to lie 80 meV below (140 meV above) the 
DP in the p (n) state. Finally, the mean free path £ s 
— v* f t is about 20 nm (40 nm) and the surface mobility 
^SdH = (e^ a )/(nfci?) is about l.OxlO 3 cm 2 /Vs (9.8xl0 2 
cm 2 /Vs) in the p (n) state. 

To infer the Dirac nature of the surface state from the 
SdH oscillations [3, 14 1, Fig. 2(e) shows the Landau- 
level (LL) fan diagram, which presents the values of 1/B 
at the rinsc-th minima in the p xx oscillations as a function 
of n osc 



4|. In the case of ideal Dirac particles, the dia- 
gram for holes (electrons) intersects the n osc -axis at —0.5 
(0.5), reflecting the it Berry phase 3; 3- However, in 
recent SdH studies of TIs ja, |B, EH exact 7r Berry 

phase has rarely been reported and this deviation has 
been a puzzle. The Zeeman coupling of the spins to the 
magnetic field can be a source of this discrepancy [Icj . 
and in addition, the deviation of E(k) from the ideal lin- 
ear dispersion also causes the Berry phase to shift from 
7r We therefore followed the scheme of Ref. 14 1 



to see if the LL fan diagram obtained for BSTS can be 
understood by considering these additional factors: The 
solid lines in Fig. 2(e) show the theoretical LLs for the 
non-ideal Dirac fermions fli| with the band parameters 
already obtained (vp = 4.6xl0 5 m/s and m* — 0.32m e ) 
and a surface g-factor g s = 20 (which is the sole fitting 
parameter). Those lines agree reasonably well with the 
experimental data for both the p and n states, supporting 
not only the Dirac nature of the observed surface state 
but also the conjecture that the holes and electrons reside 
on the same Dirac cone. 

Now we discuss the mechanism for the time evolution 
of the transport properties in BSTS. At low tempera- 
ture, the Fermi level is pinned to the IB in the bulk of 
the material [||, so the observed development of trans- 
port properties most likely comes from a change in the 
surface as schematically shown in Figs. 3(a) and (b). To 
understand the p state where holes dominate the Hall re- 
sponse, one must assume that an upward band-bending 
occurs just after cleavage, putting the Fermi level below 
the DP and simultaneously creating a hole accumulation 
layer (AL) near the surface [Fig. 3(a)]. The air expo- 
sure apparently causes n-ty pe doping on the surface as 
was reported for Bi 2 Se3 [19(, leading to a downward band 
bending [Fig. 3(b)]. 

In the above picture, there must be three transport 
channels in the p state: surface Dirac holes, bulk IB, 
and the surface AL due to the band bending. Hence, 
one may wonder if the SdH oscillations observed in the 
p state might actually be due to the AL, rather than the 
Dirac holes. Fortunately, one can see that this is not the 
case, by analyzing the non-linear B dependence of p yx . 
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FIG. 3: (Color online) (a,b) Schematic picture of the bulk and 
surface states and the surface band bending for the p and n 
states, respectively. Ec (Ey) is the energy of the conduction- 
band bottom (valence-band top), and the shaded band depicts 
the impurity band, (c) p yx (B) data measured at 1.5 K in the p 
and n states. The dashed (solid) line is the fitting of p yx (B) in 
the n-state (p-state) using the two-band (three-band) model, 
(d) Decomposition of the three-band p y x(B) fitting in the p 
state (see text). 



In the following, we discuss the analyses of the p yx (B) 
data, starting from the simpler case of the n state. 

As in BTS [5], the p yx {B) curves in the n state of BSTS 
can be well fitted with a simple two-band model described 
in Ref. jjj. The dashed line in Fig. 3(c) is a result of such 
fitting to the 1.5-K data, where the fitting parameters are 
the bulk electron density = 2.3xl0 16 cm -3 , the bulk 
mobility pt, = 190 cm 2 /Vs, and the surface mobility p s 
= 1250 cm 2 /Vs (the surface electron density was fixed at 
1.2 xlO 12 cm -2 from the SdH data). These parameters 
give the residual bulk conductivity <j& of 0.73 f2~ 1 cm~ 1 , 
and the surface contribution to the total conductance 
can be estimated as G S /(G S + cr&t) ~ 0.1, where G s ~ 
2.4xl0 -4 f2 _1 is the sheet conductance of the surface 
and t = 30 /im is the thickness. 

In the p state, the AL must also be taken into account, 
so we tried a three-band analysis in which we assumed 
that the bulk carriers are the same as in the n state. 
The solid line in Fig. 3(c) shows a result of the fitting 
to the 1.5-K data, where the fitting parameters are the 
AL mobility p s > = 770 cm 2 /Vs, the AL sheet conduc- 
tance G S ' = 2.2xl0~ 3 and the Dirac-hole mobility 
p s = 1170 cm 2 /Vs (the Dirac hole density was fixed at 
7.2xlO n cm" 2 from the SdH data). To understand the 
relative importance of the three channels, it is instruc- 
tive to examine the individual contributions to the total 
p yx : As shown in Fig. 3(d), the putative p yx {B) due 
solely to the surface Dirac holes (curve 1) is strongly 
modified when the residual bulk contribution is added 
(curve 1+2), but it is still very different from the mea- 
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FIG. 4: (Color online) Temperature dependences of p xx in 
T measured in a BSTS sample before and after reducing 
its thickness. Upper inset shows the plot of R X x(T) for the 
same set of data. Lower inset shows the high-temperature 
activation behavior in Rh(T) for the 8-pm sample, which is 
compared with the similar behavior observed in the 30-pm 
sample [Fig. 1(d) inset] shown as a dashed line. 



sured p yx (B); only when the third contribution of the AL 
is added (curve 1+2+3), the p yx (B) behavior is satisfac- 
torily reproduced. 

Based on the above analysis, one can see that it is 
impossible to interpret the SdH oscillations in the p state 
to originate from the AL: If the SdH oscillations were 
due to the AL, the third transport channel must be the 
surface Dirac holes; however, the sheet carrier density n s < 
of the third channel is estimated as n s > = G s '/(e/v) ~ 
1.8xl0 13 cm -2 , which is too large for Dirac holes for 
which kp < 5xl0 6 cm -1 and hence n s i must be < 
2xl0 12 cm -2 . Therefore, one can safely conclude that 
the SdH oscillations are due to the Dirac holes. 

Lastly, we show that the surface conductance in BSTS 
can be reasonably estimated by simply changing the 
thickness. Figure 4 shows the temperature dependences 
of Pxx of a different sample, measured first in 67-^tm 
thickness, and later cleaved down to 8-/im thick [21] . 
The overall behavior of p X x(T) is similar between the 
two curves, but a striking difference lies in their low- 
temperature saturation values pf at (for 67 /<m) and p| at 
(for 8 /jm): it is lower in the thinner sample, imply- 
ing a larger relative surface contribution. Note that the 
resistance R xx duly increases upon reducing the thick- 
ness, as shown in the upper inset of Fig. 4. Instruc- 
tively, the difference in p xx disappears at high temper- 
ature where the bulk conduction dominates. From the 
observed difference in p\ at one can estimate the surface 
and bulk contributions to the total conductivity by using 
p| at = [(G s /U) + en,]- 1 . We obtain G s w 1.8xlCT 4 ft" 1 
and <7b ~ 0.1 ri^cm" 1 , and this Ob is much smaller than 
that in the 30-/im-thick sample. This is probably because 
the number of acceptors is smaller in this second sample 
j.22| , as can be inferred in the high-temperature Rr be- 



havior [lower inset of Fig. 4]. The obtained values of G s 
and (Tb allows us to calculate the fraction of the surface 
contribution to the total conductance, G s /(G s +<Jbti), for 
the 67- and 8-pm thick samples to be 0.2 and 0.7, respec- 
tively. Therefore, when the thickness of a BSTS sample 
is reduced to < 10 pm, one can achieve a bulk TI crystal 
where the topological surface transport is dominant over 
the bulk transport. 

In summary, we show that one can achieve a 
surface-dominated transport in the new TI compound 
Bi1.5Sbo.5Te1.7Se1. 3. The surface band bending and its 
time dependence makes it possible to observe the SdH os- 
cillations of both Dirac holes and electrons, with which 
we could determine the surface band dispersion across 
the Dirac point. In addition, by analyzing the non-linear 
p yx {B) curves, we could identify the role of the surface ac- 
cumulation layer in the transport properties. Obviously, 
this material offers a well-characterized playground for 
studying the topological surface state. 
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